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Abstract 

We investigate the spectral properties of a random matrix model, which in the large 
N limit, embodies the essentials of the QCD partition function at low energy. The exact 
spectral density and its pair correlation function are derived for an arbitrary number of 
flavors and zero topological charge. Their microscopic limit provide the master formulae 
for sum rules for the inverse powers of the eigenvalues of the QCD Dirac operator as 
recently discussed by Leutwyler and Smilga. 
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The issue of chiral symmetry is fundamental in QCD. Lattice simulations indicate 
that the symmetry is spontaneously broken in the vacuum. This fact has inspired a 
large body of work in an attempt to describe the underlying mechanism. However, an 
explanation from first principles is still elusive. In a way, the spontaneous breakdown 
of chiral symmetry reflects directly on the way the quark states are delocalized in the 
vacuum |Q p], a situation reminiscent of the derealization of electrons in solids and the 
onset of conductivity [[J. 

A decade ago, Banks and Casher Q noted that the spontaneous breaking of chiral 
symmetry is related to an accumulation in the quark spectral density at zero virtuality. 
In other words, as pointed out by Leutwyler and Smilga ||, the spacing of the eigenvalues 
A is not proportional to \j\fV± as in free space, but to I/V4 (V4 is the space time volume). 
This implies the existence of a microscopic limit of the spectral density in which the 
thermodynamical limit is taken for fixed values of AV4 ||. 

In the extreme long wavelength limit, the QCD partition function can be written as an 
integral over spacetime independent Goldstone modes || . This implies strong correlations 
between the eigenvalues of the Dirac operator in the form of sum rules. These sum rules 
only involve the microscopic limit of the spectral density. 

We conjecture that the microscopic correlations between the eigenvalues of the Dirac 
operator near zero virtuality are universal, and encoded in the microscopic spectral density 
and its fluctuations. Our conjecture is motivated by the following observations: First, the 
sum rules hold for the massive Schwinger model 0. Second, they are obeyed by the Dirac 
operator in a liquid of instantons || with correlations induced by the fermion determinant 
0. Third, a random matrix model |J can be constructed that in the the large N limit 
reduces to the low energy limit of the QCD partition function for a given value of the 
vacuum angle, and thus satisfies all microscopic sum rules. Fourth, in chaotic systems 
microscopic fluctuations are universal and can be mapped on the invariant random matrix 
ensembles |8|, |9| [P!D| . Fifth, in the theory of mesoscopic systems, the Hofstadter model 
has a spectral density that coincides with that of the above random matrix model in the 



quenched approximation [fll 



The above observations suggest that a random matrix model with the general sym- 
metries of the QCD partition, might be key in understanding the general implications 



of chiral symmetry breaking in the QCD vacuum. The aim of the present paper is to 
construct the microscopic limit of its spectral density and correlations thereof. They con- 
stitute the master equations for the Leutwyler-Smilga spectral sum rules. We restrict 
ourselves to the simplest possible case of zero total topological charge. The model is out- 
lined in the next section. Using the orthogonal polynomial approach to random matrix 
theory, we derive an explicit expression for the spectral density in the chiral limit for an 
arbitrary number of flavors and zero topological charge. 

Consider a system of n = N/2 zero modes and n anti-zero modes with interaction 
given by the n x n overlap matrix T. For Nf flavors, the partition function that reflects 
the chiral structure of QCD is given by JT2|, [13], || 



Z = fvTP i T) ! n<*(% £), (1) 

where the integral is over the real and imaginary parts of the matrix elements of the 
arbitrary complex matrix T i.e., T>T is the Haar measure. In agreement with the maximum 



entropy principle [14| the distribution function of the overlap matrix elements P(T) is 
chosen Gaussian 

P(T) = exp(-^TTt). (2) 

The symplectic structure is a manifest consequence of chiral symmetry, and implies that 
the quark eigenvalues occur in pairs. The density of the zero modes, N/V4, is taken equal 
to 1, which allows us to identify the space-time volume and the number of zero modes. 

The order parameter in the study of the spontaneous breaking of chiral symmetry is 
the quark condensate defined through 

T lf =<q f q f >= lim lim -^--^—hgZ, (3) 

J J J m^OTV^oo Ndm f v ' 

where the chiral limit is to be taken after the thermodynamical limit. By writing the 
determinant as the product Yl(X^ + mj) one obtains the Banks-Casher formula [|J 

S = 7T < p c (Q) > z , (4) 

where the average < • • • >z is with respect to the partition function (1). The continuum 
spectral density is defined by 

p c (X) = lim lim -Lj(A), (5) 

all mflQ N—too jv 



and the spectral density p(A) is 

p(A)=£5(A-A n ), (6) 

where the eigenvalues ±A n are the nonzero eigenvalues of the overlap matrix in the chiral 
limit. As has been shown in || the parameter a can be identified as a = 

The sum rules recently discussed by Leutwyler and Smilga || for the QCD Dirac 
operator using chiral perturbation theory involve the average of the sums 

£ NP +q \n\m 

which can be rewritten as 



oo 



dx p -4±, rdzdyPiMpM, (8) 

o x 2 p Jo xPyi v ' 

respectively, where we have introduced the microscopic spectral density defined by 

In this letter our primary interest is to derive analytical expressions for the average of 
ps{x) and its pair correlation function that summarize the sum rules. 

The partition function (1) can be evaluated by rewriting the matrix integration in 



polar coordinates. For an arbitrary complex matrix we may write JL5 



T = UAV~\ (10) 

where U and V are unitary matrices and A is a positive definite diagonal matrix. Since 
the r.h.s has N more degrees of freedom than the l.h.s., one has to impose constraints 
on the new integration variables. This can be achieved [|15j by restricting U to the coset 
U(N)/U(1) N , where U(1) N is the diagonal subgroup of U(N). The Jacobian of this 
transformation, that depends only on the eigenvalues A& of A, is given by 

•7(A)=n(**-*«) 2 n**- (ii) 

k<l k 

The integrations over the eigenvalues and the unitary matrices decouple. The latter 
only result in an overall irrelevant constant factor and can be ignored. The eigenvalue 
distribution is thus given by 

n n 

Pn (A 1; • • • , A n ) = J(A) I] + m )) ex P("- E X D- (12) 

f k ° k=l 



The spectral density pi(A) is obtained by integration over the remaining n — 1 eigenvalues 



Pi(A] 



Yl d\ k p n (\ 



A„). 



(13) 



k=2 



These integrals can be evaluated with the help of the methods developed by Dyson, Mehta 
and Wigner (see |§ for references). The main ingredient is to write the product over the 
differences of the eigenvalues as a Vandermonde determinant, i.e. 

2 



1 



n^-A?) 2 

A 2(n-1) ... A , (n _x) 

which up to a constant can be rewritten in terms of orthogonal polynomials P k as 



fn-l(Ax) •■■ Pn-l(^n] 

The Pk will be chosen orthogonal according to the weight function 

r d(A 2 )(A 2 + m 2 )^ exp(-4A 2 )P fc (A 2 )P,(A 2 ) = 5, 
Jo a z 

For m = these polynomials are well known, 



(14) 



(15) 



kl- 



(16) 



Pk(s) 



n 



1/2 



(17) 



\a 2 F(N f + k + 1) ) " -(7* 

where the L fc ' are the associated Laguerre (Sonine) polynomials. 

The determinants can be expanded according to Cramers rule. All integrals can be 
performed immediately by orthogonality and, up to an overall constant, we are left with 

n - 1 k\ 



k% r ( N f + k + l) 



L^(z)L^(z)z N f +1 / 2 exp(-z) 



18) 



where z is defined by 



nX 2 



(19) 



The sum can be evaluated exactly with the Christoffel-Darboux formula, resulting in 



Pi (A) 



n 



nl 



a- 



T(N f + 



n. 



(L N n U{z)L N n Lt\z) - L N n f{z)L N n L + 2 \z)) z N f^eM-z), (20) 
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which, up to a normalization constant, constitutes the exact spectral density of the model 
(1). The microscopic limit is obtained by taking iV — > oo while keeping NX = x fixed 
(remember that n = N/2). This can be achieved from the asymptotic relation 

where J a is the ordinary Bessel function of degree a. The result for the microscopic 
spectral density is 

p s (x) = ^(J^(Ex) - J Nf+1 (Zx) J Nf -i{Vx)). (22) 

The normalization constant follows from the asymptotic behaviour of the Bessel function 
and the Banks-Casher relation. 

This formula reproduces all diagonal sum rules of Leutwyler and Smilga. e.g., the sum 

? a^aF (23) 

can be converted into an integral over the microscopic variable x = XN resulting in 
roo p s ( x )dx = / £ \ 2p T(2p-l)T(N f -p+l) 

Jo x 2 p \2j T{p)T{p+l)T(N f +p) { ' 

The above spectral density thus summarizes all sum rules and yields new sum rules, e.g. 
sum rules for noninteger values of p. The spectral density for 0, 1 and 2 flavors are shown 
by the full, dashed and dotted curves in Fig. 1. 

The non-diagonal sum rules require correlations in the spectral density. These cor- 
relators can be calculated using similar techniques. Also, the above arguments can be 



extended to finite quark masses [16 



The two point correlation function P2{X, A') follows from (12) by integration over n — 2 
eigenvalues and subtracting the disconnected part. Using again the properties of the 
orthogonal polynomial and the Christoffel-Darboux formulae we have 

p2{XA) -{^T(n + N F )) {ZZ) 6 { ) (25) 

with z = nX 2 /a 2 and z' = nX' 2 /a 2 . The microscopic limit p{x,y) = \im^r^ 00 p2(x/N,y/N) 
is obtained by using similar arguments to the ones used for the spectral density. The 
result is 
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p(x,y) = S xy I ^— — I (26) 



The microscopic limit ( 26| ) can be used to derive off-diagonal sum rules for the QCD Dirac 
operator. Indeed, 

< (? ^kJ > - < s m > - < ? m >! = I dxd/ W (27) 

We were not able to evaluate this integral analytically, but numerically the result is given 
by 

V 4 

(28) 



16i\^(JV/ + l)' 

for Nf = 1, 2, . . . , 20 with an accuracy of better than 1 part in 10 5 and agrees completely 
with Leutwyler and Smilga ||. 

We have shown that the microscopic spectral density following from the Gaussian chiral 
ensemble of complex matrices reproduces all microscopic sum rules for the QCD Dirac 
operator. The corresponding random matrix model reflects solely on chiral symmetry. We 
have strong indications that the microscopic spectral density is universal and should be a 
solid property of the QCD Dirac operator near zero virtuality. The overall spectral density 
is not. It would be interesting to see how the former compares to the true QCD spectral 
density following from lattice simulations. Our approach is relevant for discussing similar 
issues related to QCD in two and three dimensions, as well as the spectral properties of 
strongly coupled QED. Finally, we note that the spectral density for zero flavors describes 
the spectral correlations in the Hofstadter model as used in the framework of universal 



conductance fluctuations ITT], [17], [rj|. The physical analogy is striking, emphasizing once 
more the universality of the bulk structure of chiral symmetry breaking in the QCD 
vacuum. 
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Figure Caption 

Fig. 1. The microscopic spectral density ps(x) as a function of the microscopic variable 
x for (full line), 1 (dashed line) and 2 (dotted line) flavors. 
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